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Abstract. We are interested in the cost-convex potentials in optimal mass 
transport theory, and we show by direct and geometric arguments the equiva- 
lence between cost-subdifTerentials and ordinary subdifEcrentials of cost-convex 
functions, under the assumptions AO, Al, A2, and A3W on cost functions 
introduced by Ma, Trudinger, and Wang. The connectivity of contact sets of 
optimal transport maps follows as a direct corollary. Our approach is quite 
different from the previous result of Loeper which he obtained as the first step 
toward his Holder regularity theory of potential functions, and which was based 
upon approximation using the regularity theory of Ma, Trudinger, and Wang. 
The result in this paper improves his result, by relaxing certain geometrical as- 
sumptions on the domains and targets; it also completes his Holder regularity 
theory of potential functions on the round sphere, by making it self-contained. 



For a real valued function u on R", its subdifFerential (or subgradient) du{x) at 
X is defined locally as 

du(x) = {p G TjjR" I u{z) > u{x) + p ■ z + oQx — z|) as z ^ x}. 

It is clear that each du{xo) is a convex set and if u is differentiable at xq, du{xo) = 
{Vu(xo)}. If It is a convex function then for p £ du{xQ) the hyperplane Hp{x) = 
X ■ p -\- u{xq) supports u globally, i.e. u > Hp on R". In this paper, motivated by 
optimal mass transport theory, we are interested in the generalization of this fact 
for generalized convex functions, which are called cost-convex (or simply c-convex) 
functions, defined as the generalized Legendre transformations with respect to the 
cost function c = c{x, y) on some product space f2 x A. These are called c-Legendre 
transformations in Section [21 

Ordinary convexity is the special case c(x, y) = —x-y, and Brenier and others 
|RR| [KnSj |CuP| [Raj [Mclj showed that convex functions give potential functions for 
the optimal mass transport problem. Namely, given two mass distributions 

i.e. two Borel probability measures on R", assuming the measure /i+ does not 
charge sets of lower dimension, the unique mapping Tq which minimizes the average 
cost functional 
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among all Borcl measurable maps which push- forward /i+ onto /i_ , in the sense that 
fi-{E) — ii+{T^^{E)) for all Borel set E C R", is actually given as the gradient 
mapping of some convex function </>, i.e. Tq — V0. Here Tq and its uniqueness 
are defined in the /i+-almost everywhere sense. Brenier's result was generalized to 
other cost functions on Euclidean space [Ca2j |GMc| . to the Riemannian distance 
squared cost on manifolds [Mc2j and to the Heisenberg group [AR] , by Caffarelli, 
Gangbo and McCann, McCann, and Ambrosio and Rigot respectively. Sufficient 
conditions for the map to be smooth were investigated by Delanoe (for n = 2) 
[Dlj . CaffarelH (for n > 2) [Calj|Ca3] . and by Urbas [H], for the EucHdean distance 
squared, and then by Cordero-Erausquin and Delanoe for the Riemannian distance 
squared on the torus [Cr| and nearly flat manifolds [D2j . and finally for other cost 
functions in a series of papers by Ma, Trudinger, Wang |MTW| [TWl] . and Loeper 
[L] . We refer the reader to the book V by Villani for an introduction to optimal 
mass transportation theory and its applications in partial differential equations and 
geometry. 

For a c-convex function </> on the domain one may define the c-subdifferential 
d'^4>{xo) at xq to consist of the gradient vectors VMj^(a;o) of the functions My{x) = 
—c{x,y) + const, which globally support u at xq, meaning u{z) > My{z) for all 
z G f2 plus the equality u{xq) = My{xo). The mapping from xq to the contact 
set, namely the set of such points t/ G A, is called the c-contact mapping of and 
generalizes Brenier's subgradient mapping for ordinary convex functions. Unlike 
90, is a global notion depending on the whole domain of c, x A, but for 
convex and c(a;, y) = —x ■ y, d'^<j) = d(f> and one may naturally ask whether this 
indeed holds for other more general cost functions. 

Question 1.1. For which cost functions c on domain f2 x A, is there the equiva- 
lence between the c-subdifferential d'^cj) and the ordinary differential d(j) for c-convex 
functions 4>? 

An answer was given by Loeper [L]: there is such an equivalence d'^4> — d<j) un- 
der certain conditions on the function c, introduced by Ma, Trudinger, and Wang 
[MTWjfTWT] . who called them AO, Al, A2, and A3W , and for bounded do- 
mains ri,A C R" satisfying certain geometric restrictions with respect to c, c*, 
where c*{x,y) = c{y,x), namely that and A are mutually uniformly strictly c, 
c*-convex and c, c*-bounded (see Section [3]). Loeper's result is based upon an 
approximation argument which uses Trudinger and Wang's C'^- regularity theory of 
potential functions of optimal transport for general cost functions [TW1| , and those 
conditions on c and on SI, A which Loeper needed were actually the ones originally 
imposed by Trudinger and Wang when they applied a continuity method to the 
Monge- Ampere type equations satisfied by the potential functions. Loeper's result 
on c-convex functions was actually the first step toward his Holder regularity theory 
of potential functions for optimal transport with rough mass distributions, which 
need not be absolutely continuous with respect to the Lebesgue measure, under a 
strict version ASS of the condition A3W (see [Lj). He also extended his Holder 
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regularity theory to the standard n-dimensional sphere with cost given by its Rie- 
mannian distance squared, and to the reflector antenna problem, which helped to 
motivate Wang's interest in the theory of optimal mass transportation [Wlj |W2] 
and has also been studied by many others, e.g. [CaGHj [KoO| . 

The sufficient condition A3W of Ma, Trudinger, and Wang for regularity was 
shown by Loeper [L] to be necessary for the continuity of optimal maps with smooth 
integrable mass distributions whose Radon-Nikodym derivatives are bounded away 
from and +00. It is also necessary for the equivalence — d(j) for c-convex 
functions. 

The goal of this paper is to provide a direct and geometric proof of the equivalence 
d'^(j) = d(f> between the cost-subdifferential and the ordinary subdifferential for cost- 
convex functions cj), under Ma, Trudinger, and Wang's conditions AO, Al, A2, 
and A3W. Our result only assumes the mutual (not necessarily strict) c and c*- 
convexity of and A. For example we do not assume them to be bounded, and this 
allows us to deal with such domains as R". In fact we even assume slightly less (see 
Remark lG.ip . and our result can easily be applied to domains in the round sphere S"", 
including S"' itself (see Example 17.51 and Example 1 7. 6p . thus it completes Loeper's 
Holder regularity theory of potential functions for optimal mass transport maps on 
domains in R" or in S*" without relying on Trudinger and Wang's regularity result 
or any other method. As a consequence, it strengthens the results of Cafferelli, 
Gutierrez, and Huang [CaGH| . by weakening their assumptions on the incoming 
intensity in the reflector antenna problem. 

In this research announcement, we adopt the point of view that our domain and 
target A are either subdomains of Euclidean space, or else manifolds equipped with 
Riemannian metrics which allow us to identify their tangent and cotangent spaces. 
Unless otherwise stated, computations of higher derivatives are always carried out 
with respect to Riemannian exponential coordinates centered at the point where the 
derivatives are evaluated. In a subsequent paper, we abandon this artifice and give 
a fully covariant development of the results for general cost functions on smooth 
manifolds |KiMc| . Aside from providing conceptual clarity and geometric insight, 
this point of view offers the possibility of extending the regularity theory of Loeper 
(and presumably that of Ma, Trudinger and Wang) to other geometries. 

The present manuscript is organized as the following: In Section [2] we list the 
basic definitions and assumptions we need in this paper; in Section [3] we give the 
precise statements of main results; in Section |4] we summarize some of the rele- 
vant prior results by Ma, Trudinger, Wang, and Loeper; throughout Section [5] and 
Section [5] we prove our main results; in Section [7] we discuss some examples of 
interest. 

Acknowledgment. We thank Gregoire Loeper, Neil Trudinger, and Xu-Jia Wang 
for their helpful discussions and comments, and their recent preprints. In partic- 
ular, while writing this article, we learned from Trudinger and Wang that they 
had independently relaxed some of the restrictions on their earlier theory using a 
different method |TW2] . There they indicated the desirability of a direct approach 
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to such questions, which we beheve the present manuscript provides. We are grate- 
ful to Adrian Nachman and the 2006-07 participants of Fields Analysis Working 
Group, for the stimulating environment which they helped to create. This research 
was supported in part by Natural Sciences and Engineering Research Council of 
Canada Grant 217006-03 and United States National Science Foundation Grant 
DMS-0354729. 

2. Definitions 

In this section we list some (modified) definitions introduced by Ma, Trudinger, 
Wang, and Loeper [MTWj |TW1| [L] . In many cases, we follow the notation of [L], 
despite the fact that it is at odds with the notation of |GMc| . Let il, A be two (not 
necessarly bounded) domains in Riemannian manifolds M and N respectively; for 
example M = N = R" with the Euclidean metric. 

Definition 2.1. A cost function c : x A ^ R is said to satisfy 
(AO)ifceCit,(nxA), 

(Al) if for all x e f2 the mapping y —Vxc{x, y) is injective on A and for all 
y A the mapping x — > — Vj,c(a;, y) is injective on Q, 
(A2) if det Dlyc{x, y) 0, for all a; £ and y e A . 

Definition 2.2. (c-Exponential Maps) Under the assumption Al, we define the 
c and c*-exponential maps and their domains by 

c-Exp^ p = y iip^ -V.j;c{x,y), 

dom(c-Exp^) = {p G T^fl \ p = —Vxc{x,y) for some y G A}. 

and 

c*-Expy q = x li q = -Vyc{x, y), 

dom(c*-Exp2,) = {(7 G TyK \ q = —\Iyc[x, y) for some x € Vl}. 

The assumption A2 implies that c and c* -exponential maps are local diffeomor- 
phisms. 

Throughout this paper we will assume AO, Al, A2 for the cost function c : 
f2 X A — > R and we define c* by c*{y, x) = c{x, y). The following simple observation 
will be useful. 

Lemma 2.1. (Symmetry) For rj G T^^, £_ G TyA, 

(Dlycix, y))-^rj • ^ = 77 ' {DI^x, y)r^i 

i.e. 

[d c-^a;pJ(e-S2;pj-i(y)(77) ■i = Tl-[d c*-Sxpj^](,._^^p^)-i(,)(e) 
where '• ' means the Riemannian inner product. 



Proof. The proof is straightforward. 
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Definition 2.3. (c-sectional curvature and A3W, ASS) Define the c-sectional 
curvature for (x, y) G fi x A, 77. ^ G Txi~l, by 

(f 

at^ t=o 

for p — — Va;c(x,?/) = [c-Exp^]~^(?/). We say that c satisfies A3W (resp. ASS) if 
there exists a constant Co > (resp. > 0) such that for all {x, y) G x A, 

(2.1) &c{x,y){v,0>Co\vW for all 77 ±e 

Definition 2.4. (c-convexity) We say A is c-convex with respect to a; G H if 
(c-Exp^)~^A is convex as a set in the Euclidean space Txfl. We say A is c-convex 
with respect to if it is c-convex with respect to all x G f2. We define similarly the 
c* -convexity of Q. 

Definition 2.5. (c-convex function indexed by a set) We define the c-Legendre 
transformation L'^v of w : A ^ R U {+00} to be the function on Q given by 

L''{v){x) = sup -c{x, y) - v{y). 

J/6A 

We define a function on f2 to be a c-convex function indexed by A, if 

(2.2) (p = L'^i'v) for some function v on A 
and 

(2.3) for each x ^ there is j/ G A such that (j){x) -f v(y) — —c{x, y). 

Setting c*{x,y) = c{y,x) we may define the c*-Legendre transformation 
from functions on Q to the functions on A, and c*-convex functions indexed by 
similarly. 

Remark 2.1. ()2.3p is not explicitly assumed for c-convexity in [MTW| [TW1| [L] 
since the domains there and A are assumed to be bounded, and hence p.Sp is 
automatically satisfied due to compactness of f2 x A and continuity of c. 

Definition 2.6. (c-contact mapping) Let u be a c-convex function on fl indexed 
by A. We may continuously extend it to the closure H. of fi. We define the set- valued 
mapping G„ from to A by 

Gu{x) := {y G A I u{x) + L'^ {u){y) = —c{x, y)} for each a; G fi. 

Gu is called the c-contact mapping (or simply the contact mapping) of u, and the 
set Gu{x) is called the contact set of u at x. We can also define in a similar way 
the contact map from A to f2, for a c*-convex function w on A indexed by $7. 
Then, it is clear that for any xGil, y ^ K\i v = (u), 

G^{x) = (G„)-i(x), G„(y) = {Gu)-\y). 



Definition 2.7. (subdifFerential and c-subdifferential) Assume that VI is given 
a Riemannian inner product ■ and the distance function d induced by this metric. 
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Let M be a RU{+oo}-valued function on fl. The subdifferential du is defined locally 
for each a; G by 

du{x) = {p e Tj.Q, I u{z) > u{x) + p ■ exp^"'^ z + o(d{x, z)) as z ^ x}. 

Note that the subdifferential is always a convex set though it may possibly be 
empty. However, for a c-convex function u indexed by A, this set is nonempty, and 
we may further define the c- subdifferential d'^u for each xq d Q hy 

d^uixo) := {p e T^^il I p = -Vx=xoC{x,y), y e G'„(a;o)}. 

One should note that d'^u{x) C du{x) for any x €0. and that d'^u depends globally 
on the sets il, A. 

3. Main Results 

Using the definitions given in the previous section, we can now precisely state 
our main result. It was formulated by Loeper [L], who proved it under various 
technical hypotheses described in the section to follow using a result of Trudinger 
and Wang [TWl] . 

Theorem 3.1. (A3W CSIS) Let f2, A he two (not necessarly bounded) domains 
in Riemannian manifolds M and N, resp. (e.g. M = N = R", with Euclidean 
metric). Assume that the cost c : il x A — > R satisfies AO, Al, A2 and that A is 
c-convex with respect to f2, and SI is c* -convex with respect to A. Then the following 
are equivalent. 

1. (A3W) c satisfies A3W. 

2. (CSIS: c-subdifFerential is subdifferential) For any c-convex function (j) 
on SI indexed by A, 

d^cl) = d<j>. 

From this theorem, we can deduce the connectivity of contact sets of c-convex 
functions when we assume A3W. 

Corollary 3.2. (connectivity of contacts sets) Under the assumption of The- 
orem [Ql if further A3W holds, then for any c-convex functions (j) on indexed 
by A, the contact set G^{x) <Z K of (j) is c-convex (and so connected) for all x G f2. 

The following corollary is useful, especially when Q, — R" (see Example 17.11 and 
Example 17. 3|) . and it is parallel with Caffarelli's result |Calj for Euclidean distance 
squared cost, which required the convexity only of the target domain A. 

Corollary 3.3. (arbitrary subdomains) Put the same assumption on c, Vl, A 
as in Theorem \3.1\ . Then for any arbitrary subdomain f2i C S7 and for a subset 
Al C A c-convex with respect to f2i, and for any c-convex function (j) on ^li indexed 
by Al, we have d'^cj) — d4>. 



Proof. The proof is given at the end of Section [5l 



□ 
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4. Previous results: Ma, Trudinger, Wang, and Loeper 

We now discuss the previous results of Ma, Trudinger, Wang, and Loeper. The 
following conditions were critical in their initial studies |MTW| [TW1| [L] . Our ap- 
proach avoids many of these restrictions, so that our arguments in Section [5] and 
Section [6] only assume (not necessarily uniformly strict) c, c*-convexity — an im- 
provement which, like Corollary 13.31 was attained independently in the parallel 
work of Trudinger and Wang [TW2j , 

We denote 

(4.1) A{x,p) = —Dl^c{x,y) for y such that —Vxc{x,y) = p. 

Definition 4.1. (Uniform strict c, c*-convexity) We say that H. is uniformly 
strictly c* -convex with respect to A if its inverse c*-exponential image [c*-Expj^]~^ (SI) 
is uniformly strictly convex with respect to y S A, i.e., if is c*-convex with respect 
to A, dil G and there exists a constant 6o > such that 

[Di-fj{x) - Dp^Aijix,p)-fk]nTj > So 

for all X e dfl, p G — Va;c(a::, A) C T^Q, unit tangent vector r and outer unit normal 
vector 7. Here A is given in (14. ip . We similarly define the uniform strict c-convexity 
of A with respect to il. 

Definition 4.2. (c, c*-boundedness) We say that is c* -bounded, with respect 
to A if there exists some function (p G C^(ri), satisfying for some constant 61 > 0, 

[Dl^ix) - Dp,Aj{x,p)Dkp]C,^, > <5i|^|2 

for all X G SI and for all p G —Vxc{x,A) C T^^. Here A is given in (|4.ip . We 
similarly define the c-boundedness of A. 

Remark. In [TW1| . Trudinger and Wang mean by c-convexity and c-boundedness 
of Q, the c*-convexity and c*-boundedness in this paper. We follow the convention 
given in [L]. 

Remark. ( |TW1| . pl9) Any sufficiently small ball A is c-bounded with respect to a 
bounded set fi. 

In [MTW| [TW1| . Ma, Trudinger and Wang obtained the following regularity 
result for c-convex functions as solutions to optimal transport problems. 

Tileorem 4.1. (Ma, Trudinger, and Wang. [MTWj |TWlj see also ^ ,pl3) 
Let f2, A be two bounded domains in R". Let c satisfy A0,A1,A2, and A3W on 

S7 X A. Assume that fl, A are strictly uniformly c,c* -convex with respect to each 
other, either A is c-bounded with respect to or Q is c* -bounded with respect to 
A. Then for two probability measures po, pi on fl, A having densities pi G C^(Sl), 
P2 G C^(A) which are bounded away from 0, and for any c-convex function (j) on ft, 
with G^{fl) C A, (G^)^iiQ ~ pi, we have (f> G C^{fl). 

Using this regularity result, Loeper showed the equivalence between c-subdifferentials 
and ordinary subdifferentials (under the same assumptions of Trudinger and Wang), 
by approximating c-convex functions by smooth ones. 
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Theorem 4.2. (Xoeper [L],pi5, Theorem 3.1.) Under the same assumption on 
the cost c, domains f2, A as in Theorem \4-.l\ (especially with c or c* -houndedness) , 
we have for any c-convex function (f> on Q, with G0(ri) C A, that 

d'^4'{x) = d(j){x) for all x e fl, 

and so the set d'^(j){x) £ is convex for all x G 17. 

5. Proof of Theorem 13.11 

In this section, we illustrate a blueprint of the proof of Theorem 13. II 
We first start with the geometric setting of the proof. Pick Xm G VI and two 
different points yo,yi G A and their inverse c-exponential images po,pi G T^^^il , 
Po / Pi- Let pg — (1 — 0)po + Opi, 6* G R, be the line in T^^^fl passing through 
Pq with direction pi — Pq. For such pg G dom(c-Exp^ ), consider its c-exponential 
image yg, i.e. 

(5.1) pe = -yxc{x„i,yg) e Tx^^, ye = c-Fixp^^ {pg). 

Then, by the condition A2, we see that -^yg ^ 0. If such yg is defined at least for 
G [0, 1], then we denote by [yo, yi]xm the curve {yg}o<g<i and call it the c-segment 
with respect to Xm joining j/q to yi [L]. Such [yo)2/i]xm exists, for example if A is 
c-convex with respect to Xm . 

Definition 5.1. (Sliding and double mountains) Given Xm G Vl and j/qj Vi G A, 
define the function fg on fl, which we call the sliding mountain between yoj2/i 
centered at yg and normalized at a;„, by 

(5.2) fe{x) = -c{x,yg) + c{xjn,ye) 

where yg is the point corresponding to G [0,1] in the c-segment [yo,yi\x,n- We 
call max[/o,/i] the double mountain of j/Oi Vii normalized at x,-a. 

Definition 5.2. (The set Sg and S*^) We define for the sliding mountain /g, 

the level set Sg -.^ {x ^ Q. \ -^feix) = 0} and 

the super-level set := {x e \ -^fe{x) > 0}. 

In this notation we have the following key observation whose proof is straight- 
forward. 

Lemma 5.1. (DASM Monotonicity) For Xm G fi, suppose that yo,yi G A 
is joined by the c-segment [yo,yi]x^ with respect to Xm- Then (1) is equivalent to 
(2), and (V) is equivalent to (2'): 

(1) (DASM: Double Mountain Above Sliding Mountains) 

(5.3) /e(x) < max[/o,/i](2;), /or aZ/ (0, x) g]0, l[xa 

(2) (Monotonicity) The super-level set monotonically increases in 9, i.e. 

(5.4) 5-+ C S+ , for any < 01 < 02 < 1. 

(1') (Strict DASM) !i5.3\) holds, and the inequality is strict unless x = a;,„. 
(2') (Strict Monotonicity) holds and Sg^ n Sg.^ = {xm} if Oi < 62. 
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Remark 5.1. In this lemma, we do not assume c, c*-convexity on 51, A. 

Remark 5.2. (DASM ^ CSIS) Note that (|5.3p is equivalent to the convexity 
of the c-subdifferential 9^ max[/o, /i](a;m) at a;,„, and so to d'^ max [/o ; /i ] (^m ) — 
9max[/o, fi]{x„i)- In fact, it is easy to see as observed by Loeper (see Proposition 
2.12 in ^) that for A being c-convex with respect to 57, the property (|5.3p holds 
for every Xm,x G ft, and for all c-segments [?/0j2/i]k„, C A if and only if 9^0 = dcf) 
for all c-convex functions on indexed by A. 

Remark 5.3. (DASM ^ A3W) Loeper used Taylor expansion to prove the im- 
portant result that 

(|5.3p holds locally near Xm ^ A3W. 

(See [L] pp20-21.) 

From Lemma 15.11 Remark 15.21 and Remark 15.31 it is clear that to prove The- 
orem [331 it suffices to prove either (DASM) or (Monotonicity) for all Xm, x, 
[Uo^yilx^ while assuming A3W. We will do this by proving the following theorem 
in Section [6l 

Theorem 5.2. (A3 Monotonicity) Under the same assumption on c, il, A 
as in Theorem \3.U if c satisfies A3W then (Monotonicity) holds for all 

Xm S and 1/0,1/1 G A. Similarly A3S implies that (Strict Monotonicity) holds 
for all Xm S ^ and yo, yi e A. 

We emphasize that our proof of Theorem 15.21 thus the direction (1) =J> (2) in 
Theorem 13.11 is elementary and geometric. It is not based on the previous results 
of Ma, Trudinger, Wang, and Loeper. 

We finish this section by giving the proof of Corollarv l3.31 

Proof of Corollarv \3.3[ It is clear that the property (DASM) (|5.3p has monotonic- 
ity with respect to the domain: For fii C 51, Ai C A with A, Ai being c-convex 
with respect to 51, 51i, respectively, if (|5.3p holds for every Xm,x G 51, and for all 
c-segments [?/Oi?/i]a;„ C A, then the same holds for the pair 51i, Ai. Therefore we 
can conclude using Remark 15.21 and Theorem 13.11 that d'^(j> = d(j> for all c-convex 
functions on 51 1 indexed by Ai. 

6. Proof of Theorem 15.21 (A3 => Monotonicity) 

We divide the proof into three steps. We will first set up the necessary notation 
and show some preliminary geometric facts in the subsection 16.11 Then, we will 
show in the subsection 16.21 that the rate of expansion of 5^ along Se (see Defi- 
nition [O] and Definition 15. 2p has the same sign as the quantity -^fe on Sg. In 
the subsection 16.31 we will show the non-negativity of this expansion rate to con- 
clude the desired monotone dependence of the rising region S'^ on 9 & [0, 1] (see 
Proposition 16. 2p . and it will finish our proof of Theorem l5.2l 
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6.1. Geometric preliminaries. We first recall the definitions of pe,ye from Sec- 
tion [5l For each ye, we define 

(6.1) qe ^ ~Vyc[x„i,ye) €TygK, x„ c*-Expj^^ (ge) 

We denote ye = ^J/e, Ve — ■jffiVe^ ^^^d we also use • to denote the Riemannian 
inner product. Our calculations will be performed in geodesic normal coordinates 
of a; S O and ye G A. 

Now, let's discuss the geometric properties of Se- We first note that 

-Da:^/e •■?=(- DI^^c{x, ye) S) ■ ye, for any tangent vector ^ at x, 

and by A2 that Dxjgfe is nowhere vanishing on fi. Thus we see the zero level 
set Se of -jgf is a smooth (n — l)-dimensional submanifold in the n-dimensional 
domain $7. It is clear that Xm G Se and that 

(6.2) Se = c*-Exp^^((ge + We) n dom(c*-Exp^J), 
for the [n — l)-dimensional subspace 

(6.3) We^{w^ Ty,A \ w ■ ye ^ 0} c Ty,A. 

We parametrize by w the points in We C Ty^A and for each w let x^, be the point 
in Se <Z given by 

(6.4) w ^ -Dyc{x^,ye) + qe- 

At each point x^ G Se, denote by pe(xw) the inverse c-exponential image of ye 
given by 

(6.5) Pe{xw) = -Dxc{x^,ye), 
and denote 

(6.6) peixw) = -i7^Pe{xw)- 

aO 



Note that w - ye = 0, and by Lemma 12711 pa fa;,,, ) is orthogonal to Se at each x^j, i.e. 
(6.7) pe(a;„,) • i:>„,x,u_, = 0, 

when the derivative D^x^ is taken for any direction of w. In particular, we see 
that Se is orthogonal to pi — po at Xm for all 0. 

As our key ingredients, we consider for each 9 the functions (on Q) 

'^fs{x) = { - Dyc{x,ye) + Dyc{xm,ye)) ye, 

-J^fsi^) = ( - ^Iv^^i^-' ye) + Dlyc{x^, ye))ye ye 
+ { - Dyc{x,ye) + Dyc{x,n,ye)) ye- 
Note that at Xra, fe{xm), ^/e(x„), -^jfeixm) = 0, and for all <^ £ T^„,^, 

Dx-^fe{xrn)£. = - DxDlyc{xm,ye) C yeye 
- Dl^c{xm,ye) £, ye- 
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And, by differentiating (1 — 6)po + 9pi — —Dxc{xm,y0) twice with respect to 9, we 
see 

(6.8) = -DyyD.xc{xni, ye)yeyB - Dlyc{x„^,ye)y9 
and that 

(6.9) D,-^fe{xm) = 0. 

6.2. The rate of expansion of S'^ along So has the same sign as the 

.2 

quantity -^fe on Sg. This resuh conies form the following lemma, by letting 

a—t 

Lemma 6.1. Let g : (— e, e) xU Jl be a function, where U is an open domain 
in R". Denote gt{x) — g{t,x) and suppose the gradient S7xgt is non-vanishing on 
U for t G (— e,e). Then the moving level set Lt '.— {x £ U \ gt{x) = 0} is locally 
parametrized for small t by X{t,x), x £ Lq, where we define X{t^x) by solving the 
following first order ODE : 

[ ix(M)-- ||vf;g*(g)p V.g.(x(t,x)) 

\ X{Q,x) = X, 

where we denote gt ^9t- Moreover, the expansion rate of the super level set 
Lf := {x £ U \ gt{x) > 0} along Lt in this parametrization is given by 

Wxgtix)\\ dt Wxgt{x)\\ 

Proof. Since 

j^9t{X{t, x)) = gtiX{t, x)) + Vxgt{X{t, x)) ■ j^X{t, x) = 0, 

it is clear that this X{t,x) is indeed the desired parametrization of Lf . (|6.10p is 
clear from the construction of X{t, x). □ 



6.3. The positivity of -^fe - Now we will show that -^fi 

Se 



> under 

A3W. Under ASS condition, the equality holds only at Xm, so it is the only 
common point of the expanding boundaries Sg . This result together with the result 
in the subsection 16.21 will finish the proof of Theorem l5.2l 

Proposition 6.2. (Positivity of on Sg) Suppose that the cost function c 

satisfies AO, Al, & A2. Assume further that f2 is c* -convex with respect to A and 
that the c-segment {ye}o<0<i with respect to A. Then we have: 



(1). under A3W, 



> . 

Se 



(2). under ASS, {x) > for x ^ x^. 

Se 

To prove this proposition, let's consider the following key lemma whose proof is 
given at the end of this section. 
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Lemma 6.3. (Sliding mountain and c-curvature) We recall the c-sectional 
curvature 6c in Definition \2.3\ and the parametrization w of the set We of tangent 
vectors in Ty^A orthogonal to ye as given in (16. 3p . Let he the points in Sg C fl 
defined as in (j6.4p . and let pe, pe he the vectors defined in (|6.5[) and (j6.6p . Fix 
a tangent vector r/ g T^We, and let be the directional derivative in w of the 
direction rj and let be the second iteration of ■ Then 



Di^-j^feixw) = &c{xw,ye){D.u,x^,pe{xw})- 



Proof of Proposition \6.SX Denote 



Ve {qe + We) n dom(c*-Exp ) 



and note that Se = c*-Expy^ (Vg). We may regard the function -^fe as a 

Sg 

function, say /, on Vg by pulling it back by c*-Exp . We apply Lemma 16.31 then 



under A3W and the orthogonality (|6.7p . we see that the function / is convex. 

Now the c*-convexity of f2 at each ye implies that Vg is convex, and so / is 
a convex function on a convex domain. Then we see by (|6.9p that the function 



■jg^fe has the minimum value at Xm so implying (1), and this minimum is 

Sg 

unique under ASS when / is strictly convex so implying (2). □ 



Remark 6.1. The proof of Proposition 16.21 is the only place in the proof of Theo- 
rem [521 where we use the c*-convexity of il with respect to A. What we actually 
need here is that a dense subset of the set {qg + We)r\dom{c*-'Expy^) can be reached 
by line segments from qe in dom(c*-Expj,^). We can use this slight weaker condi- 
tion, to deal with such examples as c(x, y) = —■^ log \x — j/p in R" x R" \{x — y} 
(Example [LS]); c{x,y) = ±^\x ~ y\P,p ^ for p = -2 or -i < p < 1 (- only) 
(Example [73]) ; the cost c{x,y) = — ^logja; — ?/p restricted on S*"^^ C R" which 
appears in the reflector antenna problem (Example 17. 6p . 



In the rest of the section, we prove Lemma 16.31 



Proof of Lemma \6.3\ (Sliding mountain and c-curvature). We will use the same 
notation of Lemma 16.31 and perform our calculations in the geodesic normal coor- 
dinates of Xw and of yg. Differentiating (|6.4p twice with respect to w we see that 



(6.11) = -Dl^Dyc{xw,yg)D^Xw D^x^ - Dl^c{x^,ye)D'^^Xw 



COST-SUBDIFFERENTIAL 



13 



We compute by using (|6.8p and ()6.1ip 

2 

= -Dl^Dyyc{xw,ye)D^Xw D^x^ ye ■ ye 
- D^DyyC{x^,ye)Dl^x^u ye ■ ye 



ye 



(6.12) 



-Dl^DlyC{x^,ye)D DwXw ye ■ ye 

■ Dj;Dyyc{xw,ye) iDl^^c{xnj,ye))'^{-Dl^Dyc{xw,ye)DwXn, D^x^ 



ye ■ ye, 



To compute &^{xw,ye){Du]Xyj,Pe{xu])), we consider the line segment pe{x^ 
"0 Peixw), e R, in Tjj^ri, and its c-exponential image y(V') as 

(6.13) pe{xw) + Pe{xw) = -D.j,c{xn,,y{^))- 

Note that 2/('0)|^_q = and by differentiating (|6.13[) with respect to ip that 



Pe{xn,) = -Dlyc{xy,,ye)—y{ii^) 



therefore, we see 



ye = ^y(V') 



ip=0 



and by differentiating (|6.13p twice with respect to ip that 



= -DlyD^c{xn,,y{'4!)) ye ye - Dlyc{x^,y{^lj))—y{'ip) 



(6.14) 

Now we compute using ()6.14p 

&cixw,ye)iDn,,x^,P0{x^)) 

Dl^c{x^,y{ilj))D^Xiu Dy,x. 



^=0 



i/'=0 



DwXw ye ye 
(f 

DyDl^c{xw,ye) D^x^ D^x^ -J^vW 



(6.15) 



= -DyyDlJx^^yg) D^x^, D^Xnj ye ye 

~ DyDl^c{xw,yg) Dy,Xw D^x^ {DlyC{x^,y0)y^{-DyyDxc{x^,ye) ye ye) 

Therefore, comparing ()6.12|) and (I6.15P by interchanging the order of the deriva- 
tives and by Lemma [2T| we see that 

Diw-^feixw) = &c{xw,ye){DiuX.u],pe{xni)). 
This finishes the proof of Lemma 16.31 □ 
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7. Examples 

In the following, we demonstrate our result in some examples of cost functions 
among the ones given by Ma, Trudinger, and Wang [MTWj [TWlj . and by Loeper 
[L] . A common feature of these examples is that they only require the c-convexity of 
the target domain A with respect to the source domain fi, not the other direction, 
for the equivalence d'^(j) = d4> of c-subdifferential and ordinary subdifferential of 
c-convex functions on il indexed by A. 

Example 7.1. c{x,y) — \x ~y\'^ + \,f{x) ~ g{y)\'^ satisfies AO, Al, A2, A3W (resp. 
ASS) for il = A = R", if /, 5 : R" R convex (resp. strictly convex) and 
l^/li 1^.9! < 1- Such and A are c, c*-convex with respect to the other, since the 
domains of c, c* exponential maps are the whole space R" (One can verify this for 
example, by showing that the map —Vxcix, ■) is a proper map.) So by Corollarv l3.31 
we see that for any domains fli, Ai C R", Ai being c-convex with respect to ili, 
and for any c-convex function (p on fli indexed by Ai, we have d'^(f) = d4> without 

the c*-convexity of Qi. This result is in fact easy to show by directly computing 

,2 

-^fe and checking its positivity. 

Example 7.2. We can conclude the same as in Example 17.11 for the cost function 
c{x, y) = ^l + \x~y\^ on R" x R". 

In the following examples, we will use the notation of Section[5] These examples 
concern cost functions that have singularities. 

Example 7.3. c{x, y) = -\ log \x — on R" x R" \ {x = y} satisfies ASS. In this 
example, one can show that the curve {yejeeR is a circle (possibly with infinite 
radius) passing through x^, tangent to pi — po at x^, and that S'g's are (n — 1)- 
dimensional spheres (with different radius and centers, possibly with infinite radius) 
passing through yg and Xm, with normal vector pi — po at Xm- This geometric 
configuration shows directly the monotonicity of , as in Theorem 15.21 and thus 
either for fl — A = R" or for any C R", and A C R", c-convex with respect to 
ri, we have d'^cf) — d(j) for all c-convex functions on indexed by A. The same 
result can be obtained applying our general theory, especially Theorem 15.21 and 
Remark |6. II Note that in this case the domains dom(c-Exp^), dom(c*-Expy) (for 
= A = R") are R" \ {0} and so 

{qe + We) n dom(c*-Exp^J = [qg + Wg) \ {0}, 

(see (|6.2l)(|6.3p for the definition of We). The set {qg -f We) \ {0} is convex except 
when {0} C + Wg. In the latter case, still a dense subset of it can be reached by 
line segments from qg within it, if the dimension n > 3, as we needed in Remark l6.1l 
The case n = 1 is trivial, and the case n = 2 can be handled by an approximation 
argument. 

One may be concerned with that when poiO,pi are collinear in the tangent 
space at Xm , the corresponding c-segment is no longer differentiable when pg passes 
through the origin since the point yg in the c-segment goes to infinity. This 
problem can be bypassed by approximating the segment [po;Pi] by a line segment 
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which does not pass through the origin. Then our desired DASM survives under 
the approximation, and so it holds for il = A = R". 

Example 7.4. c{x, y) = ^^\x — y\^TP ^ satisfies A3W for p = ±2 or p = — i (— 
only) and ASS for — i < p < 1 (— only). The case p = 2 is well known and the 
same conclusion as in Example 17. II holds. For other p, the cost is not differentiable 
for X = y, but this case is similar to Example 17.31 and by the same approximation 
method using Remark l6.1l we get the same conclusion. 

Example 7.5. (The round sphere) We consider the cost c{x, y) = ^(P{x, y) on S"", 
with d the Riemannian distance function of the standard round metric of 5", i.e. 
with the diameter tt. Loeper showed [Lj that c satisfies ASS by directly computing 
the cost-sectional curvature €s- In this case, dom(c*-Expy^) = B-^[Q) C Ty^S^ 
(for n = 5") and so {qe + We) n dom(c*-Expj^J is convex (see dOldO]) for the 
definition of We) and we can apply Proposition 16.21 for those points where the 
function c{-,ye) is differentiable, to get Strict DASM on the set S*" \ Uo<e<iye, 
where we denote by ye the antipodal point of ye- Now, we use the continuity of the 
distance function to extend the Strict DASM to the whole 5". From this result 
together with Rcmark l5.2[ we can conclude that if either = A = 5" or A C 5" is 
c-convex with respect to £7 C S*" then d'^cj) = dcf) for all c-convex functions (p on 
indexed by A. 

Example 7.6. (Reflector antenna problem) The cost c{x,y) =^ — ^log|x — 
restricted on the unit sphere S*""^ C R" satisfies ASS on 5""^ x 5""^ \ {a; = y}- 
In this case, for ^l — 5""^, 

dom(c*-Exp^) = TyS""-^ R"-i 

and c*-Expy(Ty5"~^) = S"^^ \ {y}, c*-Expy(0) = y, where we denote by y the an- 
tipodal point of y in 5""^. So by our general theory (Theorem l5 . 21 and Remark l6.1|) 
as applied in Example 17.31 we can make the conclusion that if either = A = 5" 
or A C 5" is c-convex with respect to il C S", then d'^cf) = dcj) for all c-convex 
functions ^ on f2 indexed by A. 
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